CMB polarization as a direct test of Inflation 



David N. Spergel* 

Department of Astrophysical Sciences, Princeton University, Princeton, New Jersey 08544 

Matias Zaldarriaga^ 
Department of Physics, MIT, Cambridge, Massachusetts 02139 



: 

ON . 
ON ■ 

Of 

m : 
m : 

(N : 
> 

00 . 

in ■ 

o : 
r» . 
On ■ 

^: 

Oh; 

6 ■ 



X 



We study the auto-correlation function of CMB polariza- 
tion anisotropies and their cross correlation with tempera- 
ture fluctuations as probe of the causal structure of the uni- 
verse. Because polarization is generated at the last scattering 
surface, models in which fluctuations are causally produced 
on sub-horizon scales cannot generate correlations on scales 
larger then ~ 2°. Inflationary models, on the other hand, 
predict a peak in the correlation functions at these scales: its 
detection would be definitive evidence in favor of a period of 
inflation. This signal could be detected with the next gener- 
ation of satellites. 

Temperature anisotropies in the Cosmic Microwave 
Background (CMB) are one of the best probes of the 
early universe. CMB measurements are likely to answer 
one of the fundamental questions in cosmology: the ori- 
gin of the fluctuations that formed the galaxies and the 
large-scale structure. If the fluctuations are consistent 
with current models of structure formation, then accu- 
rate measurements of the fluctuations could lead to a 
precise determination of a large number of cosmological 
parameters 

There are two competing sets of theories for structure 
formation: defect models, where a symmetry breaking 
phase transition generates seeds that form sub-horizon 
scale density fluctuations, and inflationary models, where 
a period of superluminal expansion turns quantum fluc- 
tuations into super-horizon density perturbations. A fun- 
damental difference between these two mechanisms of 
structure formation is that only inflation alters the causal 
structure of the very early universe and is able to create 
correlations on super-horizon scales. 

The COBE satellite observed correlations on angles 
much larger than that subtended by the horizon at de- 
coupling (Oh ~ 1.1°) in the CMB temperature. This 
does not; however, imply that there were correlations 
on super horizon scales at decoupling because a time de- 
pendent gravitational potential will produce temperature 
fluctuations at late-times, the integrated Sachs Wolfe ef- 
fect (ISW). For example, cosmic string or texture models 
predict that most of the fluctuations observed by COBE 
were produced at z < 10. 

Measurements of temperature fluctuations at small 
scales have been suggested as a potential test of inflation: 
inflationary models and most non-inflationary ones pre- 
dict different locations and relative heights for the acous- 
tic peaks Q . Unfortunately, causality alone is insufficient 
to distinguish inflationary and non-inflationary temper- 



ature power spectra: causal sources that mimic exactly 
the inflationary pattern of peaks can be constructed ||. 
While the predicted CMB fluctuations of the current fam- 
ily of defect models differ significantly from inflationary 
predictions Q, it is useful to have model independent 
tests of the causal structure of the early universe. 

Polarization fluctuations are produced by Thomson 
scatterings during the decoupling of matter and radia- 
tion. Thus, unlike temperature fluctuations, they are 
unaffected by the ISW effect. Measurements of the polar- 
ization fluctuations are certain to probe the surface of last 
scatter. Hence, the detection of correlated polarization 
fluctuations on super-horizon scales at last scattering are 
a definitive signature of the existence of super-horizon 
scale fluctuations, one of the distinctive predictions of 
inflation.El 

We will work in the initially unperturbed synchronous 
gauge, where the metric is given by ds 2 = a 2 (r)[—dT 2 + 
(Sij + hij)dx % dx 3 }. We will consider only perturbations 
produced by scalar modes and will solve the Einstein 
equations in the presence of sources (e.g., defects) using 
the stiff approximation Q . The sources are characterized 
by their covariantly conserved stress energy tensor 0^. 
Before recombination, matter and radiation act as a very 
tightly coupled fluid, so the evolution of fluctuations can 
be described by 
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where 0/3 is the pressure, <5r and vr describe the energy 
density and velocity of the photon-baryon fluid and 5c is 
the energy density of cold dark matter. In synchronous 
gauge, the cold dark matter has zero velocity. The sum 
over N is carried out over all species and cs is the sound 



x In this letter we will consider the correlation function in 
real space (ie. as a function of the separation angle) rather 
than the usual power spectrum. By doing so we can easily 
express the causality constraint, while it would become a set 
of integral constraints that the power spectrum has to satisfy 
in the now more usual treatment in term of Czs. 
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speed. Temperature and polarization anisotropies seen 
on the sky today depend on 5f> and vr at decoupling. 

Equations ([!]) imply that the photon-baryon fluid prop- 
agates information at the speed of sound and thus can- 
not generate correlations on scales larger than the sound 
horizon. Causality on the other hand implies that the un- 
equal time correlators of the sources (6 /Jl/ (r, r)0 M „(O, r')) 
vanishes if r > r + r'. In the absence of initial 
correlations, these two conditions together imply that 
(X| r .(ni)Jf| r .(n 3 )) = if 0i2 > 29 h ~ 2°, where 
X = 5ft, vr,c^vr and r* is the conformal time of de- 
coupling. 

In the thin scattering surface approximation, equations 
are solved up to recombination and then the pho- 
tons free stream to the observer. The final temperature 
anisotropy in direction n on the sky is 

T(n) = ^| T „-n-v R | T , --J drkjii^. (2) 

The first two terms are evaluated at the last scatter- 
ing surface and the third term is an integral along the 
line of sight, the ISW effect. In non-inflationary models, 
the first two terms cannot correlate temperature fluctu- 
ations at separations larger than 29h ~ 2° but because 
anisotropies can be created later through the ISW effect 
these models can have temperature correlations on larger 
angular scales. 

Polarization is produced by Thomson scattering of ra- 
diation with a non zero quadrupole moment at the last 
scattering surface. The photons that scattered off a 
given electron came from places were the fluid had ve- 
locity vr and thus because of the tight coupling the 
photon distribution function had a dipole Ti = n • vr. 
Furthermore gradients in the velocity field across the 
mean free path of the photons (A p ) created a quadrupole 
T-2 = Xpn l n 3 diVRj in the photon distribution "seen" by 
each electron. This quadrupole generates polarization 
through Thomson scattering. 

Linear polarization is described by a 2 x 2 traceless ten- 
sor fully specified by the Q and U Stokes parameters || . 
These parameters depend on the direction of observation 
n and on the orientation of the coordinate system per- 
pendicular to n, (ei,e2) used to define them. Two inde- 
pendent combinations with spin ±2 provide a more con- 
venient description, Q±iU. Under rotation of the (ei, §2) 
basis by an angle ip, this combinations transform as 
(Q'±iU') — exp(=F2i?/>) (Q±iU), and can be expanded in 
spin ±2 harmonics, (Q±iU)(h) = J2i m a ±2,im ±2Y l m (n) 
0. An equivalent expansion using tensors on a sphere 
can be found in ||. The scattered radiation field is 
given by (Q + iU) = -3/4cr T J dfi747r(m • n')T 2 (n') cx 
Xpir^vaP d^Vj \ Tr , where ax is the Thomson scattering 
cross section and we have written the scattering matrix 
as P (m, n') = — 3/4<7T(m-n') 2 with m = ei+ie2 . In the 
last step, we integrate over all directions of the incident 
photons n'. 



As photons decouple from the baryons their mean free 
path is growing very rapidly, so a more careful analysis 
is needed to get the final polarization M , 

(Q + iU)(h) » 0.17Ar,mWait>j| T> (3) 

where At* is the width of the last scattering surface and 
is giving a measure of the distance photons can travel 
between their last two scatterings. The appearance of m 
in equation (^|) assures that (Q + iU) transforms correctly 
under rotations of (ei,e2). Equation (||) shows that the 
observed polarization only depends on the state of the 
fluid at the last scattering surface. No correlations can 
be present in the polarization for separations larger than 
~ 2° in non-inflationary models. 

Polarization can be decomposed into two distinct parts 
@: E and B with a E ,i m = ~{.a2,i m + a 2 ,i m )/2 and 
a.B,irn = —(02,1m - a 2 ,i m )/2i . Only four power spec- 
tra are needed to describe the full radiation field, three 
autocorrelations C X i = 2Tfi^ m ( a *x.im a x.im) for X = 
T,E, B and the cross correlation between E and T, 

Cci — 2T+1 Ylm( a TJm a E,lm) ■ 

Density perturbations only contribute to E We 
can illustrate this point using equation (|J) and for con- 
venience choosing n = z and m = x + iy. In the small 
scale limit, we have V 2 B = (Sg - dl)U + 2d x d y Q § 
(where V 2 = (<9 2 + <9 2 ) is the 2-d Laplacian). This gives 
W 2 B = — V 2 (z • V x vr) and is zero because the velocity 
field produced by density perturbations is irrotational. 

The correlation functions of Q and U can be defined in 
a way which makes them independent of the coordinate 
system [|| , given two directions in the sky we first rotate 
(ei, e 2 ) in each direction so that both ei are aligned with 
the great circle connecting the two points. We then use 
the Q and U as measured in this system to compute the 
correlation functions which depend only on 8 the angle 
between the two directions. They are given by, 

C (Q,U) (e) = £ aj + lpfS.*) FHe) _ C (B,B) F 2 {d)] (4) 
I 

where 

±2Y l 2 (0,^) = V '(2/ + 1)/ '4tt [F l 1 {8)±F?(0)] exp(2z0). 
Both correlation functions receive contributions from the 
E and B channels. The E channel contains all the cos- 
mological signal if there are no tensor or vector modes. 

We computed both C^' U '(9) for the model proposed 
by Turok which has a clever choice of source stress energy 
tensor that is able to reproduce the pattern of peaks of 
inflationary standard CDM (sCDM) || . The results are 
shown in figure [l]. We see that the inflationary model 
is able to produce correlations on angular scales larger 
than ~ 2°, while the other model cannot. On smaller 
angular scale than shown in figure [j], the two correlation 
functions coincide. The difference between the two mod- 
els is a result of the causal constraints and is insensitive 
to source evolution. It is also worth pointing out that 
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in inflationary models the large scale polarization is sup- 
pressed relative to the small scale signal, so we are after 
a small effect. 

Next, we estimate the expected uncertainties in mea- 
suring C^' U ^(9). Since receiver noise is the likely to be 
the dominant source of variance, we can make a sim- 
ple estimate of the total noise: it is proportional to the 
number of independent pairs of pixels, N p , at a given 
separation, 9. For an experiment with a full width at 
half maximum of 9f whm , N p = 1/2 x (4n/6j whm ) x 

(2ir9/9f w hm)- If a (Q u) i s ^ ne n °i se m the polar- 
ization measurement per resolution element, then the 
noise in the cross correlation is given by AC^ ,U ^ = 
^2jN~ p a 2 QU) ra 20 WpVO^/^y^ where 

Wp 1 = (?Iqv) Slpix/4-ir. 

We can make a more accurate determination of the 
noise using the covariance matrix of the different power 
spectra fij: 

Cov((7 ( 2 BiB)i ) = ^_(C {E , B)l + wp l e l *°lf (5) 

which give the following variances for the correlation 
functions, 

(AC (Q,U ))2 = ^ ( ^tl )2{Cov(C 2 EB)i) [F i m * 
I 

+Cov(C^ E)l ) [F 2 (9)} 2 }. (6) 

Figure |l| shows the noise in each correlation, in the 
limit where the variances are dominated by receiver noise 
(AC®) 2 = (AC U ) 2 and agree perfectly with our previ- 
ous estimates. If either the cosmic variance is important 
or the power spectra of E and B differ, then the approx- 
imate estimate of the previous paragraph is not accurate 
and the full calculation should be used to estimate the 
noise. 

The noise in the correlation functions can be reduced 
by focusing on the iS-like piece of the polarization. The 
noise m the both C( Q ' U ^(9) receives contributions from 
the variances in both the E and B spectra, but by com- 
puting both contributions separately we can show that 
the variance in E{B) makes the dominant contribution 
to AC Q (AC U ). If we filter the maps to pull out only 
the E component, then we remove not only the B signal 
but also some of the noise, AC U goes down almost by a 
factor of ~ 4. The assumption that most of the signal is 
in the E channel can be checked within the data as both 
E and B contributions can be measured separately from 
the maps. 

For the MAP satellite, without filtering the noise, 
AC( Q ' a ) w 0.36 /V9 fiK 2 , so it will not be sensitive 
enough to detect this signal, even if we combine all of 
the three highest frequencies. However, if we filter the 
map to extract the E channel signal, then the noise in 
the MAP experiment drops to AC U ~ 0.1/y/O LiK 2 , and 



the C u signal should be detectable. The PLANCK satel- 
lite, with its very sensitive bolometers, should be able to 
achieve AC^M « 0.003/v^ (J-K 2 and should easily be 
able to detect both C u and C Q . As cosmic variance is not 
the dominant contribution to the noise, an experiment 
observing a small patch of the sky could also potentially 
detect this signal. 

The temperature-polarization cross-correlation [ fL0| is 
another potential test of the origin of fluctuations: al- 
though ISW effects produce temperature fluctuations af- 
ter decoupling, we still do not expect correlations be- 
tween temperature and polarization on large angular 
scales for defect models. The correlations between tem- 
perature and polarization fluctuations directions and 
n 2 are, 

(<3(ni)T(n 2 )) = {Q\n) - \ jP drn[n{(h t3 (r)Qj) , (7) 

T* stands for the first two terms in equation (|^) and Q* 
is given by equation (||). 

In the polarization temperature cross correlation, only 
the term involving the line of sight integration could pro- 
duce correlations on large angular scales. This would 
require correlations between the late time variations of 
the metric and the velocity at last scattering. For this 
to occur in defect models, they must be moving very fast 
and remain coherent as they evolve from recombination 
to very late times. As figure 1(c) shows, even Turok's 
causal seed model, which mimics inflation remarkable 
well in the temperature correlation does not predict any 
temperature-polarization correlation. 

If gravity waves, rather than scalar modes, were the 
dominant source of the anisotropies, then they could, in 
principle, create large angular scales cross correlations. 
However, if gravity waves were significant enough to cre- 
ate a large signal, then they would be directly detectable 
in the B channel. 

Figure |l](c) shows the calculated values of the cross 
correlation together with the expected noise. The signal 
is well above the noise for MAP and Planck. The de- 
tection of a large angular scale cross correlation with no 
appreciable signal in the polarization B channel would 
put very stringent limits on the physics of models trying 
to mimic inflation. 

There is only one caveat to our argument: reioniza- 
tion. If the universe reionized very early, a significant 
fraction of the observed polarization will come from the 
rescattering of photons at late times. In most scenar- 
ios, the fraction of rescattered photons is thought to be 
less than ~ 20% pfl . Reionization has two effects on 
our argument. First, it reduces the amplitude of the cor- 
relation function by a factor exp (— 2r), where r is the 
optical depth to decoupling (r < 0.2). Second, it creates 
further structure in the correlation function on large an- 
gular scales. Fortunately, the effect of reionization can 
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be separated from that of the primordial anisotropies: it 
leaves a very specific signature in the power spectrum, 
a peak at very low I that is easily distinguish from the 
I 6 dependence expected from causality constrains alone 
@||jl4). Because of the form of F^O) and if (0), this 
peak produces an almost constant positive offset in 
and C u for angles 8 ~ 2°. Because the offset in C u is 
positive (if (0) < for ~ 2° and I < 70) reionization 
at a relatively recent epoch can never create the negative 
peak at 8 ~ 2° predicted by inflationary models. 

There is a precise signature in C^' u ' (8) on ~ 2° scales 
that would allow an unambiguous test of inflation. The 
signal is small, but within reach for the new generation 
of experiments. The cross correlation between tempera- 
ture and polarization is also expected to provide strong 
constraints that could distinguish inflation from non- ^ 
inflationary models: this signal is much larger and will , 5 
be well above the noise for MAP. The next generation "ST 
of satellites or even polarization measurements from the 
ground could provide a definitive test of the inflationary 
paradigm in the relatively near future. 
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FIG. 1. Correlation functions for Q (a) and U (b) Stokes 
parameters for sCDM and the causal seed model discussed 
in the text. The noise in their determination is shown for 
both MAP and Planck. Panel (b) shows the expected noise 
for MAP if the CMB maps are filtered to include only the 
E channel signal. Panel (c) shows the cross correlation be- 
tween temperature and polarization and the noise for MAP, 
the expected variance for Planck is even smaller. Each resolu- 
tion element in the correlation function should be considered 
independent. 
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